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Abstract 

An analytical approach has been developed to describe grand canonical equilibrium between a 
three dimensional (3D) electron system and a two dimensional (2D) one, an energy of which is 
determined self-consistently with an electron concentration. Main attention is paid to a Landau 
level (LL) pinning effect. Pinning means a fixation of the LL on a common Fermi level of the 2D 
and the 3D systems in a finite range of the magnetic field due to an electron transfer from the 
2D to the 3D system. A condition and a start of LL pinning has been found for homogeneously 
broadened LLs. The electronic transfer from the 3D to the 2D system controls an extremely sharp 
magnetic dependency of an energy of the upper filled LL at integer filling of the LLs. This can 
cause a significant increase of inhomogeneous broadening of the upper LL that was observed in 
recent local probe experiments. 

PACS numbers: 73.43. Cd 73.43.Jn 73.43.Fj 



I. INTRODUCTION 



A model of grand canonical equilibrium of a 2D system with a reservoir was proposed to 
explain the integer quantum Hall effect (IQHE) more than twenty years ago.— This so-called 
reservoir hypothesis can describe main features accompanied to the IQHE such as plateaus 
in a Hall resistance,— the Shubnikov-de Hass effect^ and oscillations of a magnetization and 
a thermoelectric power. ? Key points of the reservoir model are following: a 2D chemical 
potential is fixed by the reservoir and it doesn't depend on the magnetic field, a carriers 
concentration and 2D subband energies oscillate in the magnetic field. There was a wide 
discussion about a relevance of the hypothesis to conventional samples where the IQHE 
observed and it is still not clear. But there are structures where a reservoir is created 
nearby to a 2D system and one should take in to account the carriers exchange. These 
are tunnel and resonant tunnelling diodes.— For these devices it was shown that a carriers 
transfer between 2D and 3D systems takes place and a 2D subband energy oscillates while 
a magnetic field perpendicular to the 2D plane sweeps. 5 Moreover the 2D subband energy 
oscillations can be so strong that a pinning effect of a partial filled Landau level (PFLL) takes 
place. 6,7 Usually authors employ self-consistent calculations to get magnetic dependencies 
of a carriers concentration and a 2D subband energy and to my knowledge there is no 
any thorough analytical investigation of the LL pinning effect. In this paper I would like 
to consider thoroughly pinning of a homogeneously broadened PFLL with aims to draw 
analytical expressions, to find a PFLL pinning condition and to make other estimations. 

Recently a new local probe technique so-called a subsurface charge accumulation (SCA) 
imaging has been developed by S. H. Tessmer, R. C. Ashoori et al& to investigate a local 
compressibility of a 2D electron system (2DES). Extraordinarily sensitive to a magnetic 
field features have been revealed in the images at the magnetic fields close to the integer 
quantum Hall state of the 2DES. The SCA is measured as a charge on a tip induced by an 
applied ac voltage between the metal tip and a contact to the 2DES. A frequency of the 
signal is chosen quite low to provide the charge accumulation without delay. In this case the 
SCA is proportional to a capacitance between the tip and the 2DES that is determined by 
the local compressibility of the 2DES. In other words all applied ac voltage drops between 
the tip and the 2DES, namely, Fermi levels of the 2DES and the contact are equal. This 
means equilibrium between the 2DES and the contact. Hence at the SCA image analysis 
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one should takes in to account effects of this equilibrium. 

An analytic approach of cause requires some simplifications but the all made assumptions 
are appeared to be relevant to real tunnel structures with 2DESs. So I start in Sec. HH with 
a (5-profile of a quantum well that allows me to avoid many subbands in a 2DES. As usual 
an electron reservoir is considerably away from the 2DES that permits to neglect a variation 
of a width of a localized wave function and to use only the Poisson equation in Sec. IIHI 
instead of self-consistent calculations of the Poisson and the Shrodinger equations. The 
next assumption is in homogeneous broadening of the LLs. In this case a shape of the LL 
can be described analytically for a large LL number.— Relevance of this approach is also 
discussed in Sec. 11111 In Sec. II V I the subband energy oscillations in the magnetic field have 
been considered when PFLL pinning takes place and a starting value of the magnetic field 
has been determined for PFLL pinning. In Sec. IS I conclude this paper. 



II. MODEL OF STRUCTURE 



Let us consider a quantum well (QW) separated from a metal surface with a spacer barrier 
of thickness d (see Fig. 1). To simplify calculations a QW potential profile is supposed to 
be of the following form: 

U(z) = -a5(z). (1) 

It is well known^ that such QW has only one localized state at the energy Eq = —ma 2 /2h 2 
that is also valid in this situation because I suppose 

d > h 2 /ma. (2) 

Here 2h 2 /ma is an effective width of the localized wave function; m is an effective mass of 
the electron and H is the Planck's constant. It should be mentioned that a depletion effect 
inherent to a metal-semiconductor junction is out of the consideration and the potential 
profile of the structure at a magnetic field is shown in Figure Q As the next step I consider 
a dirty metal so that a magnetic field applied alone z axis doesn't influence on the metal 
density of states (DOS) and creates the Landau levels (LLs) in to the 2DES. This approx- 
imation is quite usual for semiconductor devices especially if one considers metal contacts. 
A temperature is supposed to be zero. 
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FIG. 1: Potential profile of the model structure with the Landau staircase and the 2D subband 
level Eq in the quantizing magnetic field. Dashed arrows show the charge transfer and the 2D 
subband shift at the magnetic field increase while the PFLL is pined at the common Fermi level. 

III. PINNING OF LANDAU LEVEL 

A magnetic field applied alone the z axis quantizes the 2D states in to the Landau levels 
(LLs). It is well known that a Fermi energy of a 2DES oscillates in a quantizing magnetic 
field for a constant electron concentration of a 2DES. In my case the 2DES Fermi level is 
fixed by the metal thus a constant concentration n is not relevant. Carried out in Ref. 
thorough calculations show that E and n are oscillating functions of the magnetic field. My 
task here is to demonstrate that there is a situation when the subband energy E changes 
in the magnetic field in such manner that the PFLL is almost fixed or pined on the common 
Fermi level of the junction. Let me start from a definition of a 2D density of states (DOS) 




H-E 



where /i is the common Fermi energy. Applying a variation 5B of the magnetic field one 
should expect the following variation of the electron concentration: 



Now I would like to substitute on the place of the 5B a variation of the cyclotron energy 
5e c = he5B/m and consider a derivative dG(e, B)/de c . To determine it's value let me 
consider the DOS in detail. This value is found as a sum of spectral functions of the 



G(e,B): 
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broadened LLs: 



(5) 



i=0 

Here (3e c is a degeneracy of the LL per unit area of the 2DES {(3 = m/nti 2 is the DOS at 
zero magnetic field). As for the spectral function it has the following dependence on e c : 



T is a width of the LL, i is a LL number, the function ip describes a shape of the LL. Here 
I neglect a LL spin splitting and an electron-electron interactions inside the 2DES. 

A shape of the LLs was a subject of an intensive study during the last three decades . 9 ' 11 ' 12 
The complexity of the problem arises from the fact that in the absence of the disorder the 
energy spectrum is discrete. As a result, the self-energy of an electron appears to be real 
in any finite order of the perturbation theory. Therefore, obtaining of a finite width of 
the LL requires summation of the entire diagram expansion. It was demonstrated 9 -*^ that 
such a summation is possible when the number of the LL is large. The simplifications, 
arising in this limit, are different in the case of a short-range and a smooth disorder. In 
the former case the correlation radius R c of the disorder less than the magnetic length I 
and only a subsequence of diagrams without self-intersections contributes to the self-energy, 
or, in other words, the self-consistent Born approximation 9 becomes asymptotically exact. 
The shape of the LL in this case is close to semielliptical. For a smooth disorder, with 
£ < i? c < i? L ^ i l / 2 l (R L is the Larmour radius) the LL shape had been found Gaussian. 
In this case the random potential can be renormalized to effective one with R c > Rl and 
hence it becomes a long-range potential.— In the case of the long-range disorder R c 3> Rl 
the semiclassical approach is valid when the disorder simply modulates the subband energy 
Eq and LL staircase follows this modulation. This means an inhomogeneous broadening of 
the LL and it will be considered somewhere else. Thus only the short-range disorder leads 
to the homogeneous LL broadening. 

As for the model under consideration the correlation radius R c could not be longer than 
the distance between the 2DES and the metal d because of the metal screening^. Hence the 
short-range disorder approximation is adequate when £ > d. In this case the LL spectral 
function can be written as 
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where F = ^ed 2 ', 7 = {2h/Tir) 1 / 2 and r is a relaxation time at B — 0. 9 

Other mechanisms of LL broadening such as phonon-electron, electron-electron interac- 
tions are negligible because the temperature is zero. LL broadening caused by electron 
tunnelling between the 2DES and the metal is possible but also negligible due to a very low 
tunnel transparency in accordance with Eq. (j2J. 

Thus having the LL shape determined, one can make the next step. From Eqs. (@J) and 
(J3J) one can get: 

Sn = —Gp SE + £ E ° (P <Pi(?, £ c) + Pe c d(p f £ ' £c) ^j de Se c . (8) 

Here I have Gf = G(fi — Eq, B). In accordance with equations (jSJ) and (j3J) one can simplify 
the equation in to the following: 

Sn = —Gp SE + -Se c + T ^ V /fe c ^£l^ de Se c . (9) 
e c Jo i=0 oe c 

According to Eq. (J2J) the derivative in the last term of this equation can be expressed as: 
dLpi(e,e c ) _ <pi(e,e e ) (% + 1/2 ^ e \ d(pj(e,e c ) ^ 



(11) 



de c 2e c \ 2 2e c J de 

or 

difi{e,e c ) _ (i + 1/2\ dipi(e,e c ) 1 d[eipi(e,e c 



de c \ 2 J de 2e c de 

Now one can easy perform integration in the last term of Eq. (jOJ) and get the following: 

6n = -G F 6E + -8e c -f2f3eJ^^ + l -±^ ^(// - E , e e ) Se c , (12) 

£c i=0 \ £ c J 

where I have already used that r C e c and set ipi(0,e c ) = 0. Following this approximation 
further one can get from Eq. (JHJ) that Gp = (3e c ipN(n — Eq, e c ) (N is a number of the PFLL). 
In this case the Eq. fll2j) is modified to the following: 

^ r-m n , ft*-E Q iV+l/2\ . , . 

5n = -G F 5E + -5e c - I ^ ^ u + -^\G F 5e c . (13) 

The next step is to find a relation between 6n and SEq. The rigorous way is to solve 
self-consistently the Poisson and Schrodinger equations but I have inequality (J2J) that allows 
me to neglect the variations of the wave function and QW widths in compare with that of 
the QW depth. The depth variation can be easy evaluated from the variation of the electric 
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field between the metal and the 2DES. In this case the subband level Eq follows the QW 
depth. Hence one can get the following: 

5E = e 2 5n/C . (14) 

Cq is a geometric capacity per unit area, i.e., Cq = K/d, where k is a permittivity of the 
structure material. Combining Eqs. (JT3jl and ()14|) one can determine the derivative dE /de c 
as: 

dE = 2e 2 n - \p - E + (N + l/2)e c ]G F e 2 

de c 2e 2 e c G F + 2e c C ' [ } 

According to this equation the derivative of a central energy of the PFLL dE^/dec can be 
calculated as 

dE ^ _ N 1 dE o _ 2e2n + 2 ^C (N + 1/2) + [{N + l/2)e c - /x + E }e 2 G F 
de c 2 de c 2e c G F + 2e c C ' [ ' 

PFLL pinning takes place when the derivative is zero or negligibly small, i.e., dE^/dec <^ 1. 
Hence the pinning condition is the following: 

G F 3> n/e c + NCo/e 2 , (17) 

here I have used that in the case of the large G F the Fermi level is very close to the LL 
center and the following inequation is justified \i — E — (N + l/2)e c < T < e c . Also the 
iV 3> 1 is supposed. 



IV. RESULTS AND DISCUSSION 

Since the PFLL is pined on the Fermi level the subband energy Eq should have a magnetic 
dependence. From Eqs. (I15|) . (|T7|) one can get a quite expectable result for this dependence: 



dE /de c = -(N+ 1/2). (18) 

This means that under the pining condition a basement of the LL staircase Eq shifts down 
on the energy scale and this shift compensates the PFLL energy increase due to the increase 
of the LL staircase step with the magnetic field (see Fig. 1). For integer filling of the LLs 
G F is zero and from Eqs. ()15j) and (|lfi|) one can get: dEo/de c = e 2 n/(e c Co) and dE\^jde c = 
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iV + 1/2 + e 2 n/ (e c Cq). In the case of integer filling factor the electron concentration can be 
expressed as n = N(3e c then one can get 

dE /de c = e 2 (3N/C , (19) 

and 

dE LL /de c = (1 + e 2 f3/C )N - 1/2. (20) 

Here I suppose the N-th LL as just emptied. Let me make evaluations of the derivatives 
for the typical structures parameters. For the AlxGai_xAs/GaAs heterostructures I have 
(3 = 0.28 x 10 11 (meVcm 2 )" 1 and d = 100 nm. Hence I get e 2 (3/C Q « 46. This means that 
the LL shift becomes 47 times faster in the magnetic field at the integer filling factor. As for 
the subband energy Eq it also has the large derivative that signifies a rapid increase with the 
magnetic field. Thus the central energy demonstrates the step-like behavior while Eq 
has a saw-tooth oscillation on the magnetic field (see Fig. EJ). This is qualitatively consistent 
with experimental observations and numerical calculations.— 

The strong field dependencies of the levels Eq and can explain an extraordinarily field 
sensitivity of the SCA features.— Actually let me imagine an inhomogeneous variation of 
the subband and the PFLL energies in the 2DES. This means different local Fermi densities 
of states Gp that corresponds to different derivatives dE^/de c and dE /de c in accordance 
with Eqs. (|15|) and ()16j) . Due to the difference of the derivatives PFLL broadening or a 
dispersion of a E^ distribution changes with a magnetic field variation. The grater are 
values of the derivatives dE^/de c the grater is variation of the dispersion. Thus at low 
values of Gp, namely, near integer LL's filling the E^ dispersion should be most sensitive 
to the magnetic field variation. Moreover since the derivative dE^/de c increases at the Cq 
decrease, the closer is the reservoir to the 2DES (the lesser is d) the lesser are the derivatives 
dEu^/dec [see Eq. (JTHj) ] and the lesser is the sensitivity of the Gf and the E^ dispersions 
and the SCA image contrast. This was observed in SCA images measured on a 2DES tunnel 
coupled to a n + reservoir J 4 For instance, one can compare Figure 3 in Ref. Q and Figure 3 



in Ref. ll3P^. In former case the image contrast changed markedly at the magnetic field 
variation of 5Bi = 1 kOe and in the last case it is of 5B 2 = 50 Oe. According to Eq. (J20|) 
at low values of Gp dE\^lde c is inverse proportional to Cq that can be estimated from 

n 

parameters of investigated structures. Thus in the case of Ref. |l4j the tunnel distance is 
40 nm that corresponds to Cq = n/d = 2.8 x 10~ 3 F/m 2 . The parameters of Ref. 
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FIG. 2: Magnetic dependencies of the subband energy Eq and the central energy £ll of the iV-th 
PFLL. 

thoroughly considered in Ref. where a C str ay can be considered as a capacity between the 
2DES and the ohmic contact or the reservoir. Authors had also evaluated the stray capacity 
as 1 fF. At this case a charging area of the 2DES is a disk of L = 90 nm diameter under the 
tip. Hence the equivalent capacity Cq s can be estimated as Cq s = C stTa ,y/L 2 = 0.12 F/m 2 . 
By this means one can compare the relation of the magnetic field variations 5Bi/5B 2 = 20 
and that of the LL centers derivatives, which can be estimated as a relation of the capacities 
Cos/Cq ~ 43 in accordance with Eq. ([20)1 . One can see that the relations are of the same 
order and the model proposed in this paper can be applied for description of the local probe 
experiments. However to describe the SCA features in details some improvement of the 
model is necessary. First of all the in-plane charge transfer should be considered. In this 
case the C dependence upon Gp should be found. Many-body effects also can be included 
in the model. The derivative dEp/dn therewith should be locally found and substituted in 
the place of Gp. 
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Analytic Eqs. (|16|) and f)15|) allow to determine a magnetic field value when PFLL pinning 
starts. This question is not simple because according to Eq. (fTBj) the derivative dE^/de c 
is always positive. Hence there is no a sharp LL pinning transition but some criterion can 
be found. To do this let me consider a discrepancy between a LL pinning picture and a 
LL behavior in a closed 2DES with a constant electron number. In the last case a subband 
energy Eq is constant and dE^/dSc = N + 1/2. Under the LL pinning condition E^ 
is constant and Eq. ()18j) is valid. Thus it seems to be reasonable to choose dE\^jde c = 
—dE /de c = (N + 1/2) /2 as a starting point for LL pinning. Hence according to Eq. ()15|) 
and one can get the following: 



where s c q is the cyclotron energy at the LL pinning start. Since I'm interesting in the starting 



Eq.©. Also taking in to account that N + 1/2 = v/2 — 1 w n/ f3e C Q for the half populated 
PFLL the start cyclotron energy e c0 can be found as 



It worth to note that this value is one order larger than cyclotron energy of the LLs resolution, 
which is determined from e cl w V as e c % ~ 7 2 . 

To finish the section it will be pertinent to discuss real structure parameters and magnetic 
fields to those the developed model is applicable. First of all one should satisfy N ^> 1 
and e c ^> T conditions. This definitely requires quite high mobility and concentration of 
electrons in a 2DES. Let me suppose a ten times larger than one as a much more than one. 
So for a quite usual mobility of the 2DES as 77 = 100m 2 /(Vxs) one can resolve LLs at 
B = = 0.01 T. Since according to Eq. (}2*2*|) PFLL pinning will start at B = B & 0.1 T. 
In this case one should provide at least ten LLs populated that corresponds to E F » 2 meV 
or n fh 6 x 10 10 cm -2 . It is worth noting the larger is the electron concentration the better 
is the model applicability. Now one should provide homogeneous LL broadening that is to 
locate a metal at the distance d <^il from the 2DES plane. At B = B the magnetic length 
i is about 250 nm. Thus d = 25 nm will be sufficient to assure homogeneous broadening of 
the LLs up to B = 10 T. At the last step one need to justify inequality (J2J). This means that 
QW width should be about 3 nm. 



e 2 n + e c0 C {N + 1/2) = N+l/2 
£co£ 2 Gf + £ c qCo 2 



(21) 



point I should consider Gf value in the LL center, i.e., Gp 



1 /2 

2/3e c o /7T7 in accordance with 




(22) 
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V. CONCLUSIONS 



In conclusion it is noteworthy that some results can be expanded on a closed 2DES with 
a constant electron concentration n. In this case it would be more relevant to say about a 
Fermi level capture by the PFLL. Condition of the capture can be derived in similar way 
like that of LL pinning. It can be obtained from Eq. (jl7j) setting Co = that is equivalent 
to an infinitely far reservoir. The same setting should be apply to get a start of the Fermi 
level capture from Eq. (J22J)- 

Summarizing thermodynamic equilibrium between a 2DES and a 3D electron system in 
a quantizing magnetic field has been considered in the case of the homogeneous broadened 
Landau levels. Parameters of the proposed model have allowed to consider analytically such 
effects as PFLL pinning and the subband energy oscillations in the 2DES. A condition and 
a starting point of PFLL pinning have been defined and a starting value of a cyclotron 
energy has been determined. Moreover found dependencies of and E up on Cq allow to 
describe qualitatively the extraordinary sensitivity of the SCA features observed recently in 
local probe experiments. All results are in a good qualitative agreement with experimental 
data. 
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